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Abstract 



We present an updated analysis of astro-physical solutions, two-flavor 
MSW solutions, and vacuum oscillation solutions to the solar neutrino 
anomaly. The recent results of each of the five solar neutrino experiments 
are incorporated, including both the zenith angle (day-night) and spectral 
information from the Kamiokande experiment, and the preliminary Super- 
Kamiokande results. New theoretical developments include the use of the 
most recent Bahcall-Pinsonneault flux predictions (and uncertainties) and 
density and production profiles, the radiative corrections to the neutrino- 
electron scattering cross section, and new constraints on the Ga absorption 
cross section inferred from the gallium source experiments. From a model 
independent analysis, arbitrary astrophysical solutions are excluded at more 
that 98% C.L. even if one ignores any one of the three classes of experiment, 
relaxes the luminosity constraint, or allows more suppression of the 7 Be than 
8 B flux. The data is well described by large and small mixing angle two-flavor 
MSW conversions, MSW conversions into a sterile neutrino with small mix- 
ing, or vacuum oscillations. We also present MSW fits for nonstandard solar 
models parameterized by an arbitrary solar core temperature or arbitrary 8 B 
flux. 
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I. INTRODUCTION 



The solar neutrino problem currently provides one of the most compelling experimental 
signatures for the physics beyond the standard model. The Mikhyev-Smirnov-Wolfenstein 
(MSW) effect [l| via neutrino mass and mixing provides a complete explanation of the 
existing solar neutrino data, while astrophysical solutions, even those with drastic alter- 
ations of the standard solar model, simply fail. The difficulty with astrophysical explana- 
tions persists even if we ignore data of any one of the three kinds of experiments, i.e., the 
Homestake chlorine experiment , the water Cerenkov experiments of Kamiokande P and 
Super-Kamiokande ||], or the gallium experiments of SAGE || and GALLEX ||. (The 
experimental results are summarized in Table H along with the SSM predictions fTj.) The 
successful results of gallium source experiments |||5| and the excellent agreement between 
the standard solar model predictions and the recent helioseismology data further reinforce 
our confidence in neutrino oscillation solutions ||. The new generation of solar neutrino 
experiments, such as Super-Kamiokande and Sudbury Neutrino Observatory (SNO), will 
provide critical tests of the MSW predictions for the neutrino energy spectrum distortion, 
the day-night rate asymmetry, and the charged-to-neutral current ratio. 

In this paper we examine the current status of the solar neutrino problem for astrophys- 
ical solutions, MSW solutions, and vacuum oscillation solutions. The data as of February 
1997, including the final Kamiokande results and the preliminary Super-Kamiokande data, 
are used. This is the first MSW analysis using the entire data of the Kamiokande spectrum 
and day-night asymmetry. We also incorporate the latest standard solar model with dif- 
fusion effects @, the radiative corrections for the neutrino-electron scattering cross section 
@, the improved determination of the 8 B decay spectrum ||10|| , and the constraint on the 



gallium cross section from the source experiments [11]. The calculations in this paper are 



described in detail in our previous works [|12|-[T6[|, including the model-independent analysis 
for astrophysical solutions, MSW calculations, the day- night effect, and consistent treatment 
of solar model uncertainties. We consider only two-flavor oscillations because of their sim- 
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plicity and viability. We referred to Ref. |17j for a recent analysis for three-flavor oscillations 
and Ref. ]Tj| for recent developments in neutrino physics. 

This paper is organized as follows. In Section [TT] we reexamine the general astrophysical 
solutions and show their failure with much stronger statistical significance than before. This 
is true even if we ignore any one of the three types of experiment or the solar luminosity 
constraint. We also discuss Monte Carlo evaluations of goodness of fit when the number of 
degrees of freedom (DOF) effectively becomes zero or negative. In Section |T| the constraints 
on the MSW parameters are updated. The Kamiokande spectrum result by itself excludes 
the adiabatic (horizontal) branch almost entirely. The MSW solutions with nonstandard 
core temperatures, and with nonstandard 8 B flux, and oscillations to sterile neutrinos are 



also examined. Vacuum oscillation solutions are discussed in Section |IV|. We show that 
the Kamiokande spectrum data considerably restricts the allowed parameter space. The 
conclusions of our analysis are given in Section [V]. 



II. ASTROPHYSICAL SOLUTIONS 



The incompatibility of astrophysical solutions and the solar neutrino data has been inves- 



tigated in many ways. These include the failure of explicit nonstandard solar models [19,20 



the comparison of the Homestake and Kamiokande results [21], lower core temperature fits 
|2"2| , [T2|| , and so on. One can generalize the argument against astrophysical solutions by a 
model independent analysis using pp, 7 Be, 8 B, and CNO fluxes as free parameters under 
minimal assumptions on the solar luminosity, the beta spectrum shape, and the detector 
cross sections. The details of our analysis is described in [0,0 (similar analyses are found 
in |23|-[26|1 ) . We will display the results of the fits in the 0(Be) —0(B) plane, both normalized 
to the SSM values (0(Be) S sM = 5.15 x 10 9 and 0(B) S sm = 6.62 x 10 6 in units of cm" 2 s _1 

0)- 

The constraints from individual data are shown in Fig. |I[ The combined result from 
Kamiokande and Super-Kamiokande determines the 8 B flux only. The CNO flux as well as 
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the 7 Be and 8 B fluxes are used as free parameters in fitting the Homestake result. In fitting 
to the combined Ga result from SAGE and GALLEX, the pp flux is also varied as a free 
parameter subject to the luminosity constraint. 

A comparison of Fig. [I] with our original analysis in 1993 (Fig. 1 in Ref. [I3"|) displays a 



dramatic improvement in statistics, especially in the water Cerenkov data and the gallium 
data. The addition of the high-statistics Super-Kamiokande result with about 1000 events 
from the first 100 days [] has reduced the uncertainty in the 8 B flux measurement in half. 
The low rate and the precision of the gallium result alone impose serious problems for 
astrophysical solutions. The 8 B flux allowed by the Homestake and gallium data each is 
smaller than the Kamiokande-Super-Kamiokande measurement for almost the entire range 
of the 7 Be flux. In addition the Homestake and gallium together are incompatible, since 
for 0(Be) ~ and 0(B) ~ 0, the gallium data requires the CNO flux to be ~ 0, while the 
Homestake data requires the CNO flux 4.9 times larger than the SSM value. 

The severity of the problem with astrophysical solutions can be seen by applying the joint 
analysis to all the data, shown in Fig. |^. We allow the 7 Be flux to be negative. The best 
fit of the combined observations is in the non-physical region: 0(Be)/0(Be)ssM — —0.6 ±0.4 
and 0(B)/0(B) SSM = 0.4 ±0.05 with Xmm = 0-5 for 3 data points, one luminosity constraint, 
and four free parameters: i.e., zero degrees of freedom (DOF). Within physical parameter 
space (0(Be) > and 0(B) > 0), the best fit is 0(Be)/0(Be) SS M < 0.1 and 0(B)/0(B) SS m = 
0.38 ±0.05 with Xmm = 9-2. The usual prescription for goodness of fit (GOF) evaluations f28fl 
does not apply since we have zero DOF. (Later we will encounter fits with negative DOF). 
In addition the probability distribution is non-Gaussian due to the physical constraints, i.e., 
the fluxes should be non-negative. Generalization of GOF by employing the Monte Carlo 
method is necessary. 



1 The previous Kamiokande experiment collected a total of only 600 events in 5.7 years. 

2 We can also allow the CNO flux to be negative. In this case the allowed region is essentially 
unbounded for 0(Be) <C 0. 
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GOF in this case is denned by the probability to obtain Xmm as large as 9.2 or larger by 
chance due to the experimental uncertainties, if the best fit fluxes are true. Our Monte Carlo 
construction is (1) to take the central flux values of the fit, (2) calculate the solar neutrino 
rates for the three experiments, (3) generate Monte Carlo distributions for each experi- 
ment assuming the actual experimental uncertainties (7.7, 7.4, and 9.7% for the Homestake, 
Kamiokande/Super-Kamiokande, and gallium experiments, respectively), (4) for each Monte 
Carlo data set, apply our model independent analysis to obtain the x 2 minimum. Note that 
when one has N parameters and M constraints with Gaussian errors and N > M, this pro- 
cedure can be done analytically, reproducing the usual \ 2 distribution for N — M dimensions 

The Monte Carlo distribution of the x 2 minima is shown in Fig. Q the Xmm = 9-2 from 
the actual data is also indicated. The probability of getting x 2 minimum larger than 9.2 by 
chance is 0.6%. That is, our model independent analysis excludes the best fit astrophysical 
solution at the 99.4% C.L. Q 

Next we consider the same analysis but ignoring one of the constraints. Fig. f|, ^j, and |6| 
show the results each without the Homestake data, the water Cerenkov data, and the gallium 
data, respectively. Fig. is the result with all experiments, but without the luminosity 
constraint. (Violations of the luminosity constraint would be possible if the properties of 
the solar core were somehow varying on a time scale short compared with 10 4 years.) The 
corresponding GOF of the best fit in the physical region is 98.9, 98.3, 98.9, and 98.3% C.L., 
respectively. Although the constraints on the fluxes are somewhat relaxed by ignoring one 
of the data or the luminosity constraint, the essential problem with the poor fit remains. In 



3 We have also considered two alternative nonstandard evaluations of GOF. If one assumes (some- 
what arbitrarily) that the CNO flux is fixed to zero J13|| , the fit is for 1 DOF and Xmm = 9-2 
corresponds to 99.8% C.L. If one defines GOF by the ratio of the volume of the likelihood function 
integrated within the physical region to the volume integrated in the entire parameter space (in- 
cluding negative fluxes), the ratio is 0.2%, or the exclusion is 99.8% C.L. Thus, the both estimates 
are similar to the Monte Carlo result. 
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addition one can also see the persistent problem of the strong suppression of the 7 Be flux, 
which is difficult to obtain by astrophysical effects in general. 

III. MSW SOLUTIONS 
A. Kamiokande spectrum and day-night data 

The Kamiokande experiment has completed its measurements and published the results 
of total rate, spectrum data, and day- night data [[|. The MSW parameters can be con- 
strained by those results. Note that the spectrum distortions and day-night time-dependence 
are not expected with standard neutrino physics and are powerful indicators of physics be- 
yond the standard model. The spectrum shape measured in Kamiokande is consistent with 
the one expected from the undistorted 8 B /3-decay spectrum albeit the uncertainties are 
large. The shape is inconsistent with the strong suppressions at large energies expected 
in the MSW adiabatic branch (Am 2 ~ 10" 4 eV 2 and sin 2 29 ~ fCT 4 - O.f), and the data 
exclude the region almost entirely as shown in Fig. [8[ This exclusion is simply due to lack of 
distortions in the data and independent of the uncertainties in the initial 8 B flux. The data, 
however, do not constrain the nonadiabatic (diagonal) branch, in which spectrum distortions 
are smaller. 

The day-night result was published as one day-time rate and five bins for night-time 
bins. The binning was cos# = — 0.2, 0.2 — 0.4, ... 0.8 — 1, where 9 is the angle between the 
direction to the Sun and the nadir at the detector. Within the experimental uncertainties 
the six bins are consistent and thus exclude a large region in which day-night asymmetries 
due to the Earth effect are expected. The excluded region is shown in Fig. || 

The allowed parameter space from the total Kamiokande rate is shown in Fig. || This 
constraint is model dependent and we have assumed the Bahcall-Pinsonneault model (1995) 
including its uncertainties. In this and other fits the correlation in the theoretical uncer- 
tainties between the flux components and between the experiments are included. 

Unfortunately we cannot combine the spectrum and day-night data since the errors in 
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those are strongly correlated and the correlation matrix is unpublished. Fig. |10| and |TT] show 
the allowed region when the total rate is combined with the spectrum data and day-night 
data, respectively. 



B. Preliminary results from Super-Kamiokande 

Recently the Super-Kamiokande collaboration reported the results of about 1000 events 
from the first 100 days of data M. The total rate (see Table |) is consistent with the 
previous Kamiokande rate, and new uncertainties are much smaller. When combined with 
the Kamiokande total rate, the error is reduced almost by half. The new constraint is shown 
in Fig. O The uncertainties in the MSW parameter space are now dominated by the 8 B flux 



error in the solar model calculations (~ 15%). Future measurements of model-independent 
quantities such as the spectrum shape and day-night effect and also the charged-to-neutral 
current ratio in SNO, are essential to confirm the MSW interpretation and to improve the 
determination of the MSW parameters. 

C. MSW combined results 

We next consider the MSW constraint including the Homestake and gallium data. The 
two separate allowed regions are shown in Fig. [13|. The fit includes Kamiokande day-night 
data and the averaged Super-Kamiokande data. (The allowed regions are essentially identical 
even if the Kamiokande spectrum data are used.) Both allowed regions provide a good fit. 
The x 2 minimum for 7 degrees of freedom is 5.9 (55% C.L.) and 6.4 (49% C.L.) for the 
small-angle and large-angle solution, respectively. (Details are listed in Table [II].) The 
fit for the large-angle solution improved from previous analyses (see for example Ref [jl5| ) 
due to the larger 8 B flux in the new SSM and a new, smaller Kamiokande and Super- 
Kamiokande rate, both of which reduces the relative difference between the Homestake rate 
and the Kamiokande/ Super-Kamiokande rate and allows energy-independent flux reduction 
as expected in the large-angle region. 
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We have also considered oscillations to sterile neutrinos |30]]. The GOF for the large 
angle region is 94% C.L. However, the 95% allowed region defined by x 2 < Xmm + 6-0 does 
not appear in the sin 2 29 - Am 2 parameter space (Fig. [TJ] and Table [TIT]). The large angle 



solution for sterile neutrinos is also severely constrained by big bang nucleosynthesis pO|,|3T 



While nonstandard solar models, as discussed in Section [Q|, cannot solve the solar neu- 



trino problem, the MSW effect can be also considered with nonstandard solar models ||12| ! [15 
Many of those models may be parameterize by nonstandard core temperature (Tc) or sim- 
ply a nonstandard 8 B flux, whose uncertainties might be larger than the SSM estimate. We 
consider joint fits of Tc or 8 B flux, in addition to the MSW parameters. 

When T c is used as a free parameter, the neutrino fluxes can be scaled according to the 
power law. From the Monte Carlo investigation of the SSM, the indices of the power law 



are obtained in Ref. |32|, based on the Monte Carlo SSMs by Bahcall and Ulrich [ 19| . [] The 



combined Homestake, gallium, Kamiokande, and Super-Kamiokande data constrain 

T C /T S C SM = 0.9918;°* (la), (1) 

and T C /T C SM = 0.95 - 1.02 for 95% C.L., where T^ SM is the the SSM value (1.567 x 10 7 K). 
This result is in excellent agreement with the SSM range 1 ± 0.006 \W[\: the data are 



consistent with the SSM prediction in the presence of the MSW effect. Our likelihood for 



Tc is shown in Fig. 15. The corresponding MSW parameter space is shown in Fig. 16 



Next the 8 B flux is used as a free parameter, the combined data determines 

0(B)MB) SSM = °- 76 ^3o (I*)- (2) 

and 0.31 - 1.50 for 95% C.L. Although the uncertainty is large, the result is consistent with 
the SSM range llo!i7 0- O ur likelihood for 0(B) and the corresponding MSW regions are 



shown in Fig. O and [TR 



The Monte Carlo estimate for the model with diffusion is not yet available. 
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IV. VACUUM OSCILLATION SOLUTIONS 



The simple two-flavor vacuum oscillations are still a phenomenologically viable solution 
|. Those solutions require tuning of Am 2 and the Sun-Earth distance at the 5% level to 
explain the observations, which is a conceptual setback. 

Some parameter space for the vacuum oscillation for Am 2 ~ 10~ 10 eV 2 predicts a rel- 
atively strong energy dependence, and the recent Kamiokande spectrum data alone can 
exclude a wide range of parameters, as shown in Fig. [TP]. When combined with the results of 



Homestake, gallium, and Super-Kamiokande, we find three separate allowed regions within 
a narrow range of parameters [Am 2 = (5 — 8) x 10~ n eV 2 and sin 2 29 = 0.65 — 1] as shown 



in Fig. The GOF for the best fit parameters is 9.9 for 9 DOF, which is acceptable. 
Details of the fits are listed in Table [IV] For comparison we show five allowed regions when 
the Kamiokande spectrum data are ignored Fig. |2l|) . 

V. CONCLUSIONS 

Although the general scope of the solar neutrino problem has not chanced since the first 
result of the gallium experiment in 1992, the improved accuracy of the solar neutrino data 
provide a more robust assessment of solutions. The astrophysical solutions in general have 
difficulties unless all experiments are wrong, or at least two out of three data and the SSM 
are wrong. 

The MSW effect provides viable solutions: the small mixing-angle solution (sin 2 29 ~ 
0.008 and 5 x 10 -6 eV 2 ) and the large-angle solution (sin 2 29 ~ 0.6 and 1.6 x 10~ 5 eV 2 ), 
assuming the latest SSM by Bahcall and Pinsonneault. Oscillations to sterile neutrinos are 
possible, but only for small angles. When the core temperature or the 8 B flux is used as 
a free parameter, the joint data determines those at the 3% and 30% level, respectively. 
Those ranges are consistent with the SSM predictions. Vacuum oscillations are still viable 
for Am 2 ~6x lO" 11 eV 2 and sin 2 29 ~ 0.9. 

The Kamiokande day-night data and spectrum data each exclude a large parameter space 
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for MSW, independent of SSM predictions. We expect Super-Kamiokande will provide 
those with much improved accuracy and eventually, along with the SNO neutral current 
measurement, single out the solution of the solar neutrino problem. 
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TABLES 



TABLE I. The standard solar model predictions of Bahcall and Pinsonneault (BP SSM) |7j 
and the results of the solar neutrino experiments. 





BP SSM 


Experiments 


Homestake 


9.3 t\l SNU 


2.55 ± 0.14 ± 0.14 SNU (0.273 ± 0.021 BP SSM) 


Kamiokande 

Super-Kamiokande 

Combined 


6.62 ±5;?1 a 


2.80 ± 0.19 ± 0.33 a (0.423 ± 0.058 BP SSM) 
2.51 toil ± 0.18 a (0.379 ± 0.034 BP SSM) 
2.586 ± 0.195 a (0.391 ± 0.029 BP SSM) 


SAGE 

GALLEX 

Combined 


137 t 8 7 SNU 


69 ± 10 t\ SNU (0.504 ± 0.089 BP SSM) 
69.7 ± 6.7 tH SNU (0.509 ± 0.059 BP SSM) 
69.5 ± 6.7 SNU (0.507 ± 0.049 BP SSM) 


a In units of 10 6 cm~ 2 sec~ 


-l 




TABLE II. The best fit parameters, the x 2 minimum, and confidence levels of GOF for the 
combined MSW fits. 


Small Angle Large Angle 


sin 2 26 
Am 2 (eV 2 ) 




8.2 x 10~ 3 0.63 
5.1 x 10" 6 1.6 x 1(T 5 


X 2 (7 d.f.) 
P(%) 




5.9 6.3 
45 49 


TABLE III. The best fit parameters, 
oscillations to sterile neutrinos. 


X 2 minimum, and GOF for the combined MSW fits for 


Small Angle Large Angle 


sin 2 26 
Am 2 (eV 2 ) 




1.0 x 10~ 2 0.72 
4.0 x 10~ 6 8.9 x 1(T 6 


X 2 (7 d.f.) 
P(%) 




6.7 13.7 

54 94 
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TABLE IV. The best fit parameters, \ 2 minimum, and GOF for the combined vacuum oscil- 
lation fits including the Kamiokande spectrum data. 





Solution 1 


Solution 2 


Solution 3 


sin 2 29 


0.83 


0.90 


1.0 


Am 2 (eV 2 ) 


7.9 x KT 11 


6.6 x 10 -11 


5.2 x 10~ n 


X 2 (9 d.f.) 


11.4 


9.9 


11.9 


P(%) 


75 


64 


78 
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FIG. 1. The constraints on the 7 Be and 8 B fluxes at 90% C.L. from the Homestake result 
(below the dotted line), the combined Kamiokande and Super-Kamiokande results (between the 
dot-dashed lines), and the combined SAGE and GALLEX results (below the dashed line). The 
SSM range is also shown (solid line, 90% C.L.) 
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FIG. 2. The constraints from the combined CI, Ga, and Gerenkov experiments at 90, 95, and 
Yo C.L. Also shown are the Bahcall-Pinsonneault SSM region at 90% C.L. 0, the core temper- 
ature power law and standard and nonstandard solar models including the recent 3 He diffusion 



model by Cunning and Haxton (see Ref. [16] for references for the other models). A smaller 
1S17 cross section moves the solar model predictions to a smaller 8 B flux as indicted by the arrow. 
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FIG. 3. The Monte-Carlo distribution of \ 2 minima when the best fit fluxes in the physical 
region (0(Be)/^(Be)ssM = and 0(B)/^)(B)ssm = 0.35) are assumed. The actual \ 2 minimum 
of the combined observations are indicated by the arrow. The best fit astrophysical solution is 
excluded at the 99.4% C.L. 
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FIG. 4. The result with the water Cerenkov and gallium data only. 
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FIG. 5. The result with the Homestake and gallium data only. 
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FIG. 6. The result with the Homestake and water Cerenkov data only. 
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FIG. 7. The result with the combined data but without imposing the luminosity 
constraint. 
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FIG. 8. The MSW parameter space excluded by the Kamiokande spectrum data and day-night 
data. 
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FIG. 10. The MSW parameter space allowed by the Kamiokande total rate and spectrum data. 
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FIG. 11. The MSW parameter space allowed by the Kamiokande total rate and day-night data. 
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FIG. 12. The comparison of the MSW parameter space allowed by the Kamiokande total rate 
(shaded region) and the combined Kamiokande and preliminary Super-Kamiokande rate (dashed 
lines). The theory error (~ 15%) is the leading uncertainty in the combined fits. 
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FIG. 13. The result of the MSW parameter space (shaded regions) allowed by the combined 
observations at 95% C.L. assuming the Bahcall-Pinsonneault SSM with He diffusion. The con- 
straints from Homestake, combined Kamiokande and Super-Kamiokande, and combined SAGE 
and GALLEX are shown by the dot-dashed, solid, and dashed lines, respectively. Also shown are 
the regions excluded by the Kamiokande spectrum and day-night data (dotted lines). 
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FIG. 14. Same as Fig. 13 except that this is for oscillations to sterile neutrinos. 
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FIG. 15. The likelihood distribution of the core temperature from the simultaneous MSW fit 
to the combined observations. 
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FIG. 16. The MSW parameter space allowed by the combined observations when the core tem- 
perature is used as a free parameter. The model independent exclusion regions by the Kamiokande 
spectrum and day-night data are also shown. 
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FIG. 17. The likelihood distribution of the core temperature from the simultaneous MSW fit 
to the combined observations. 
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FIG. 18. The MSW parameter space allowed by the combined observations when the 8 B flux 
is used as a free parameter. The model independent exclusion regions by the Kamiokande spectrum 
and day-night data are also shown. 
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FIG. 19. The vacuum oscillation parameter space excluded by the Kamiokande spectrum data. 
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FIG. 20. The vacuum oscillation parameter space allowed by the combined observations 
including the Kamiokande spectrum data. 
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FIG. 21. The vacuum oscillation parameter space allowed by the combined observations but 
without the Kamiokande spectrum data. 
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